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ABSTRACT 

We searched for steady PeV gamma-ray emission from the Monogem ring 
region with the Tibet air shower array from 1997 February to 2004 October 
. No evidence for statistically significant gamma-ray signals was found in a 
region 111° < R.A. < 114°, 12?5 < decl. < 15?5 in the Monogem ring where 
the MAKET-ANI experiment recently claimed a positive detection of PeV high- 
energy cosmic radiation, although our flux sensitivity is approximately 10 times 
better than MAKET-ANI’s. We set the most stringent integral flux upper limit 
at a 99% conhdence level of 4.0 x 10“^^ cm“^ s“^ sr“^ above 1 PeV on diffuse 
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gamma rays extended in the 3° x 3° region. 

Subject headings: (stars:) supernovae: individual (Monogem Ring) — gamma 
rays : observations 


1. Introduction 

In a recent observation at PeV energies, the MAKET-ANI air shower experiment at 
Mount Aragats (E44°10', N40°30'N; 3200 m above sea level) claimed a detection of signihcant 
excess (Oa) of cosmic-ray events within a 3° x 3° search window (111° < a < 114°, 12?5 < 
S < 15?5) in the Monogem ring region using the air shower data recorded from 1997 to 2003 
(Chilingarian, Martirosian, & Gharagyozyan 2003) Naturally, the signihcant excess may be 
attributed to PeV gamma rays, because the Larmor radius i?L ~ 0.4 / Z pc at 10^^ eV in the 
galactic magnetic held of 3 /iG is too small to reach the Earth without dehection compared 
with the distance of 300 pc between the Earth and the Monogem ring, and the mean decay 
length of a neutron A ~ 10 pc at 10^® eV is also too short. 

The Monogem ring is a dihuse (extended with a diameter of 25° in the sky) supernova 
remnant (SNR) associated with the radio pulsar PSR B0656-I-14 at a distance of approxi¬ 
mately 300 pc (Thorsett et al. 2003). This SNR is a bright source in the soft X-ray region. 
According to a detailed observation of the Monogem ring made by the ROSAT X-ray survey 
(Plucinsky et al. 1996), the average temperature of thermal emission is 6.15 in log(T/l K). 
If the Monogem ring is modeled as an SNR in the adiabatic stage of evolution, the initial 
ambient density is estimated to be 5.2 x 10“^ cm“^, the initial explosion energy is estimated 
to be 1.9 X 10^° erg, the radius to the shock front is estimated to be 66.5 pc, and the age is 
estimated to be 86 thousand years. The Monogem ring has been considered to be a possible 
candidate for the particle acceleration site of cosmic electrons (Kobayashi et al. 2004) and 
nuclei (Thorsett et al. 2003) in our Galaxy. 

The cosmic-ray flux we observe at the Earth roughly shows a power-law spectrum at 
energies from 10® to 10^® eV. There is a general consensus that the stochastic shock ac¬ 
celeration at SNRs in our Galaxy could explain the cosmic-ray energy spectrum up to 
approximately 10^® eV. The maximum energy is estimated to be only ^ Z x 10^^ eV by 
parallel shock acceleration where the normal axis of shock front is parallel to the direction 
of the interstellar magnetic held. On the other hand, the oblique shocks accelerate parti¬ 
cles more efficiently than the parallel shocks (Jokipii 1987). The maximum energy based on 
the oblique shock acceleration increases by several orders of magnitude, and could explain 
the cosmic-ray energy spectrum up to approximately 10^'^ eV (Kobayakawa et al. 2002). In 
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this case, the cosmic rays beyond 10^'^ eV are assumed to be of extragalactic sources, such 
as active galactic nuclei and gamma-ray burst . Among recent models, the “single source 
(SS) model” (Berezhko et ah 1996; Erlykin & Wolfendale 2004) is interesting in that a sin¬ 
gle nearby SNR, such as the Monogem ring, mainly contributes to the cosmic-ray intensity 
observed at the Earth. To explain the shape and intensity of the cosmic-ray energy spec¬ 
trum using the SS model, the most likely parameters of the single SNR should be 300-350 pc 
distant and 90-100 thousand years old (Erlykin & Wolfendale 2003). If such a strong cosmic- 
ray accelerator lies near the Earth, we may observe an anisotropy of cosmic rays from its 
direction, or we may detect high-energy photons that are emitted from these high-energy 
charged particles by the nonthermal processes. 

In this paper, we report on the search for diffuse/pointlike PeV gamma-ray emission 
based on the data recorded from 1997 to 2004 around the Monogem ring region by a large 
air shower array with a total area of 36,900 m^ constructed in Tibet. 


2. Experiment 

The Tibet air shower experiment has been successfully operated at Yangbajing (E90°31', 
N30°06'; 4300 m above sea level) in Tibet, China since 1990. The Tibet I array was con¬ 
structed in 1990 (Amenomori et al. 1992) and it was gradually expanded by 1994 to the Ti¬ 
bet II that consisted of 185 fast-timing (FT) scintillation counters placed on a 15 m square 
grid covering 36,900 m^ and 36 density (D) scintillation counters around the FT-counter 
array. Each counter has a plastic scintillator plate (BICRON BC-408A) of 0.5 m^ in area 
and 3 cm in thickness. A 0.5 cm thick lead plate is put on the top of each counter in order to 
increase the counter sensitivity by converting gamma rays into electron-positron pairs in an 
electromagnetic shower (Bloomer, Linsley, & Watson 1988; Amenomori et al. 1990). All the 
FT counters are equipped with a fast-timing photomultiplier tube (FT-PMT; Hamamastu 
HI 161) measuring up to 15 particles, and 52 out of 185 FT counters, which are arrayed at 30 
m lattice intervals, are also equipped with a wide dynamic range PMT (D-PMT; Hamamatsu 
H3178) measuring up to 500 particles. The time and charge information of the FT-PMTs is 
recorded, while only the charge information of D-PMT is recorded. All the D counters are 
also equipped with FT-PMT and D-PMT, where only charge information of both PMTs is 
recorded. An event trigger signal is issued when any fourfold coincidence occurs in the FT 
counters recording more than 0.6 particles. The mode of energy of the triggered events in 
Tibet H is ~10 TeV. 

From 1996 to 2003, we upgraded the array, and, at present, it consists of 761 FT 
counters covering 50,400 m^ and 28 D counters around them. In the inner 36,900 m^, FT 
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counters are deployed at 7.5 m lattice intervals. Since 1999 October, we have called this 
upgraded array Tibet III. The mode of energy of the triggered events in Tibet III is 3 
TeV. Using Tibet III, we observed the multi-TeV gamma-ray flares from Mrk 421 in 2000 
and 2001 (Amenomori et ah 2003a) and the cosmic-ray anisotropy (Amenomori et al. 2004; 
Amenomori et al. 2005a), surveyed the northern TeV gamma-ray sky (Amenomori et al. 2005b), 
and measured the all-particle energy spectrum (Amenomori et al. 2003b) and chemical com¬ 
position of the cosmic rays in the knee region (Amenomori et al. 2003c). 


3. Analysis 

In the present paper, we employ the data obtained by the 185 FT counters and the 36 
D counters corresponding to the Tibet II array configuration for the whole period in order 
to simplify the analysis. We collected 1.6 x 10® air shower events during 1717 detector live 
days from 1997 February 15 to 2004 October 10 after the quality cut and the event selection 
based on the following simple criteria: 

1. Air shower corelocation .—Among the three hottest counters in each event, two should 
be contained in the inner 36,900 m^. 

2. Shower size .—Spo which is the sum of the number of particles per m^ counted by the 36 
D counters and 52 out of 185 FT counters that have a D-PMT, should be more than 100. 

3. Zenithangle .— The zenith angle of the arrival direction should be less than 40°. 

To examine the performance of the Tibet II array, we use a Monte Carlo (MC) simu¬ 
lation. We employ the CORSIKA version 6.200 code (Heck et al. 1998) for the generation 
of air shower events and the EPICS uv8.00 code (K.Kasahara 2005)^ for the detection of 
shower particles with scintillation counters, respectively. The primary gamma rays are sam¬ 
pled from the Monogem ring orbit assuming a differential power-law spectrum with spectral 
index —2.0 from 10 TeV to 30 PeV. 

According to the MC simulation including the quality cut and the event selection, the 
mode energy of gamma rays is 150 TeV, the angular resolution is less than 0?3, and the 
effective area for gamma rays is nearly 2.5 x 10^ m^ for a 3° x 3° search window in a 
diffuse source analysis and 1.6 x 10^ m^ for a 0?5 x 0?5 search window in a pointlike source 
analysis. The gamma-ray energy is estimated from Spo by the MC simulation and the energy 
resolution is less than 30% above 150 TeV. The systematic pointing error is estimated to be 


^Additional information on Kasahara (2005) is available at littp://eweb.n.kanagawa- 
u. ac. j p/ ~ kasahara/ResearcliHome/EPICSHome /. 
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less than 0?02 by the Moon’s shadow in cosmic rays (Amenomori et ah 2003a). The source 
positional uncertainty is typically less than 0?1 at 150 TeV for gamma rays, assuming a 5 a 
signihcance level. The center of the Monogem ring (5 = 14°) stays in the held of view (with 
a zenith angle of <40°) for 380 days out of 1717 live days. 

Subsequently, we use the right ascension scan method to search for PeV gamma-ray 
sources that follows the same analysis method and parameters employed by the MAKET-ANI 
experiment (Chilingarian, Martirosian, & Gharagyozyan 2003). First, each event is sorted 
by its arrival right ascension and declination into a Aa x AS = 3° x 3° rectangular cell. Off- 
source events are taken from all the cells (except the on-source cell) in the same declination 
band as the on-source cell. The signihcance of the source in each cell is calculated based on 
exsquation (17) of Li & Ma (1983). 

We scanned the celestial sky in the declination band from —5?5 to 66?5 in the whole 
right ascension range 0° - 360°. We also scanned the whole Monogem ring region and the 
region around it in the declination band from 0° to 30° in the right ascension range 80° - 
130° with a Aa x AS = 0?5 x 0?5 search window analysis for a pointlike source. 


4. Results 

Figure 1 shows the number of events in each of the 120 cells in the declination band 12?5- 
15?5 by the Aa x AS = 3° x 3° search window analysis with Spn > 1000 (corresponding to > 
1 PeV). The shaded histogram denotes our actual result, and the dashed one is the expected 
excess from the MAKET-ANI result. The MAKET-ANI experiment detected a signihcant 
excess in the direction 111° < a < 114°, 12?5 < S < 15?5 at the 6 a statistical signihcance. 
But no signihcant signal was detected by the Tibet air shower array(—O.Ga). It should be 
noted that the Tibet air shower array accumulated about 100 times as much statistics as the 
MAKET-ANI experiment, corresponding to an approximately 10 times better sensitivity, 
which is easily calculated by the number of background events as summarized in Table 1. 
In addition, the signihcance distribution as shown in Figure 2 implies the absence of such a 
bright and dihuse source above 1 PeV, even in the all sky survey (0° < a < 360°, —5?5 < 
S < 66?5). 

Taking the pointing errors of both experiments into account, we scanned the area of a ± 5°, 
(5 ± 5° around the source direction claimed by the MAKET-ANI experiment by sliding the 
center of the 3° x 3° search window by 0?5 steps in right ascension and declination; however, 
no signihcant excess was found again. 

Table 1 demonstrates the energy dependence of the result in the region suggested by the 
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MAKET-ANI experiment by the 3° x 3° search window analysis, also conhrming no signal 
detection. As the MAKET-ANI experiment detected a signihcant excess in the various 
energy thresholds of >800 TeV, >1 PeV, and >2 PeV, we would have detected a signihcant 
signal at 50 ± 10 a (estimated from the number of on-source and off-source events based 
on equation (17) of Li & Ma (1983), which is a different dehnition from MAKET-ANPs) in 
at least one energy threshold even if the relative energy scale uncertainty between the two 
experiments differed by a factor of 2. 

Figure 3 shows the signihcance map of the whole Monogem ring region and the region 
around it (80° < a < 130°, 0° < 5 < 30°) based on a hner window search of 0?5 x 0?5 cells 
at energies > 1 PeV. Again, no signihcant signal was found. There are two directions with a 
signihcance > Aa at (a, 5) = (108?75, 18?75), (120?25, 12?75) that are very far from the point 
where the MAKET-ANI experiment claimed the detection of a signal. The expected number 
of directions above 4a in a normal Gaussian distribution with 6000 trials is 0.19, and the 
probability of getting more than 2 in a Poisson distribution with mean value of 0.19 is 1.6 x 
10“^. Therefore, the deviation may be due to statistical huctuations. The results of the 0?5 
X 0?5 window searches did not show any signihcant energy dependence, and no signihcant 
deviation of signihcance distribution from a normal Gaussian distribution was found at the 
other energy thresholds (> 500 TeV, > 800 TeV, and > 2 PeV). 


5. Discussions 

No evidence for statistically signihcant gamma-ray signals was found in a region 111° < 
a < 114°, 12?5 < S < 15?5 in the Monogem Ring where the MAKET-ANI experiment 
recently claimed a positive detection of PeV high energy cosmic radiation, although our hux 
sensitivity is approximately 10 times better than MAKET-ANPs. We set the most stringent 
integral hux upper limit at a 99% conhdence level of 1.1 x 10“^^ cm“^s“^/(2.66 x 10“^sr) 
= 4.0 X 10“^^ cm“^ s“^ sr“^ on steady dihuse gamma rays > 1 PeV extended within the 
rectangular region in the Monogem Ring assuming a diherential spectral index —2.0. One 
of the potential possibilities for explaining the discrepancy between the two experiments 
could be transient emission that occurred at occasions when we stopped data acquisition for 
annual maintenance, calibration, upgrading jobs, etc. Another could be the strong transient 
emission of PeV gamma rays that incidentally or periodically occurred during the 3 hr when 
only the MAKET-ANI experiment could observe because the two experimental sites are 
separated by about 45° in longitude. Third, the excess events that they detected might be 
due to statistical huctuation (Ghilingarian et ah 2005). 

No signihcant signal was found in the whole Monogem Ring region based on a 0?5 x 0?5 



window search for a pointlike source at energies > 1 PeV. We also set 99% confidence-level 
flux upper limits of 2.6 and 5.4 x 10“^^ cm“^ s“^ on the steady gamma rays > 1 PeV from 
PSR B0656-I-14 and Geminga, respectively, assuming pointlike sources with a differential 
energy spectral index —2.0. The KASCADE group also reported that no signihcant signal 
was seen in the sub-PeV region at the suggested location by the MAKET-ANI experiment 
and PSR B0656-I-14 by a pointlike source analysis (Antoni et al. 2004). Furthermore, we re¬ 
ported on the result of wide sky survey for steady TeV gamma-ray pointlike sources elsewhere 
(Amenomori et al. 2005b), although no signihcant point source was found in the Monogem 
Ring region above a few TeV and above 10 TeV. 


This work is supported in part by Grants-in-Aid for Scientihc Research on Priority Area 
(712) (MEXT) and also for Scientihc Research (JSPS) in Japan, and by the Committee of 
the Natural Science Foundation and by the Chinese Academy of Sciences in China. 
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Fig. 1.— Number of events above 1 PeV in each of the 120 cells in the declination band of 
12?5-15?5 observed by the Tibet air shower array during 1717 live days from 1997 Febru¬ 
ary 15 to 2004 October 10. The shaded histogram represents our actual result, while the 
dashed histogram denotes the number of events expected from the MAKET-ANI result 
(Chilingarian, Martirosian, & Gharagyozyan 2003). 
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Fig. 2.— Significance distribution of 2880 Aa x A6 = 3°x3° cells in the declination band - 
5?5 < S < 66?5 above 1 PeV. Each cell is independent of one another. The solid line indicates 
a normal Gaussian fit to the data. 


































































12 



85 90 95 100 105 110 115 120 125 130 


r-^ 30 

do 

^ 25 
CO 

20 


15 


10 


Geminga 


PSR B0656+14 


Entries 6000 


>1PeV 


a [deg.] 



Fig. 3.— Top: Significance map in the Monogem Ring region and the region around it above 
1 PeV. Each cell is independent of one another. PSR B0656+14 is located at the center of 
the region. The rectangular region indicates the area where the MAKET-ANI experiment 
claimed a positive detection of a signal. Bottom: Significance distribution in the region 
shown in the map. The solid line indicates a normal Gaussian £t to the data. The two cells 
with > 4 a are located at {a, S) = (108?75, 18?75) and (120?25, 12?75), respectively. 
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Table 1. Enegy dependence of the number of events in the cell 111° < a < 114°, 
12?5 < S < 15?5 and the background 


Primary Gamma-ray 
Energy*^ [TeV] 

prBG b 
^''Tibet 

juON 

i\tBG b 

MAKET 

ATON 

MAKET 

>500 

5664.7 

5605 

58 

84 

>800 

2609.6 

2618 

26 

57 

>1000 {N^ > 10®) 

1785.9 

1759 

18 

43 

>2000 

506.3 

490 

4 
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^The primary gamma-ray energy of the MAKET-ANI experiment was estimated by the 
air shower size reported in (Chilingarian, Martirosian, & Gharagyozyan 2003), where >10® 
at MAKET-ANI altitude corresponds to >1 PeV (Erlykin & Wolfendale 2004) and where we 
scaled the other energies, while that of the Tibet air shower array is estimated by the air 
shower size Spn, described in the text, in the Tibet air shower array. 

^Mean value of the background cells in the declination band 12?5 < S < 15?5. Subscripts 
“Tibet” and “MAKET” represent our result and the MAKET-ANI result, respectively. 







